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We provide a brief introduction to the high density effective theory of QCD. As an 
application, we consider the instanton correction to the perturbatively generated 
gap in the color superconducting phase. We show that the instanton correction 



^ ■ becomes large for ~ 1.25 GeV in Nf = 2 QCD, and for /i ~ 750 MeV in A^^ = 3 

_ QCD with a massive strange quark. We also study some other numerical issues 

related to the magnitude of the gap. We find, in particular, that a renormalization 
^ • group improved gap equation does not give results that are substantially different 

I ' from a gap equation with a fixed coupling. 

^ - 

1. Introduction 



■ Over the last several years we have seen rapid progress in the theoretical 

, study of very dense hadronic matter. Many new phases of strongly inter- 

im ' acting matter, such as color superconducting quark matter and color-flavor 

locked matter have been predicted i>2,3,4,5^ Reviews of these developments 
can be found in 6,7,8,9, lo^ Exotic phases of matter at high baryon density 
may be realized in nature in the cores of neutron stars. In order to study 
this possibility quantitatively we would like to develop a systematic frame- 
work that will allow us to determine the exact nature of the phase diagram 
as a function of the density, temperature, the quark masses, and the lep- 
ton chemical potentials, and to compute the low energy properties of these 
phases. In this contribution we give a brief review of an attempt to use 
effective field theory methods to address this problem. 

2. High Density Effective Theory 

At high baryon density the relevant degrees of freedom are particle and hole 
excitations which move with the Fermi velocity v. Since the momentum 
p ~ w/i is large, typical soft scatterings cannot change the momentum by 
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very much. An effective field theory of particles and holes in QCD is given 

by "'12 

L = J2^l{iv . - \g-.GI^ + ..., (1) 

V 

where — (!,«)■ The field describes particles and holes with momenta 
p = liv + l, where I <S /z. We will write I = + l±_ with /|| = v{l ■ v) and 
Z_L = Z — In order to take into account the entire Fermi surface we have 
to cover the Fermi surface with patches labeled by the local Fermi velocity. 
The number of these patches is n„ ~ {fj? / K\) where A_l <C is the cutoff 
on the transverse momenta l±_ . 

Higher order terms are suppressed by powers of 1/ /i. As usual we have 
to consider all possible terms allowed by the symmetries of the underlying 
theory. At 0{l/fi) we have 

—i^lDlA - 9^l^^^^^v^ ■ (2) 

At higher order in there is an infinite tower of operators of the form 
with w = and n = 2ni +712-1. At 0(1/^2) 

the effective theory contains four-fermion operators 

L = ^ ^ ^ C^^' {Vi ■ V2,V-i- V3,V2- V3) 

^ Vi r,r' 

• (^,,r^,,) {^PlT'^iyiv, +V2-V3- Vi). (3) 

There are two types of operators that are compatible with the restriction 
Vi + V2 = V3 + V4. The first possibility is that both the incoming and 
outgoing fermion momenta are back-to-back. This corresponds to the BCS 
interaction 

^ = ^EE^'''^r(^^-^")(^''r^-)(^I'rVl.,), (4) 

^ v,v' r,r' 

where v ■ v' = cos 6 is the scattering angle and B^^ {v ■ v') is a set of 
orthogonal polynomials. The second possibility is that the final momenta 
are equal to the initial momenta up to a rotation around the axis defined 
by the sum of the incoming momenta. The relevant four-fermion operator 
is 

^ = 7^2 E E ('^)^r' ■ ^) (V'^rv-^^) (^trvl,) , (5) 

v,v',ii>r,r' 
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where v, v' are the vectors obtained from v, v' by a rotation around Vtot = 
V + v' hy the angle (f). 

The four-fermion operators in the effective theory can be determined by 
matching moments of quark-quark scattering amplitudes between QCD and 
the effective theory. The matching conditions involve on-shcU scattering 
amplitudes in BCS and forward kinematics. The scattering amplitude in 
the effective theory contains almost coUinear gluon exchanges which do not 
change the velocity label of the quarks as well as four-fermion operators 
which correspond to scattering involving different patches on the Fermi 
surface. There are several contributions in the microscopic theory that are 
absorbed into four-fermion operators in the effective theory. Examples are 
large angle scatterings and non-perturbative instanton effects 

3. Power Counting 

In this section we briefly discuss the power counting in the high density 
effective theory. We will denote the small scale by I <C /x. We first discuss 
the scaling properties of a generic operator. We assume that v ■ D scales 
as I, tpy scales as l^^"^, scales as I, and D±^,v ■ D ^ I. Complication 
arise because not all loop diagrams scale as l^. In fermion loops sums 
over patches and integrals over transverse momenta can combine to give 
integrals that are proportional to the surface area of the Fermi sphere, 

2tt^J (27r)2 2tt^ J An' ^ ' 

V 

These loop integrals scale as t^, not Z*. In the following we will refer to 
loops that scale as P as "hard loops" and loops that scale as l'^ as "soft 
loops" . In order to take this distinction into account we define and 
to be the number of soft and hard vertices of scaling dimension k. A vertex 
is called soft if it contains no fermion lines. In order to determine the I 
counting of a general diagram in the effective theory we remove all gluon 
lines from the graph, see Fig. 1. We denote the number of connected pieces 
of the remaining graph by Nc- Using Euler identities for both the initial 
and the reduced graph we find that the diagram scales as with 

5 = Y,[{k- 4)lf + {k-2- h)V^] +Eq+A- 2Nc. (7) 

i 

Here, fk denotes the number of fermion fields in a hard vertex, and Eq is 
the number of external quark lines. We observe that in general the scaling 
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Figure 1. Counting hard loops in the effective field theory. If all (soft) gluon lines are 
removed the remaining fermionic loops contain sums over the velocity index. 

dimension S increases with the number of higher order vertices, but there 
are two important exceptions. 

First we observe that the number of disconnected fermion loops, Nc, 
reduces the power 6. Each disconnected loop contains at least one power 
of the coupling constant, g, for every soft vertex. As a result, fermion loop 
insertions in gluon n-point functions spoil the power counting if the gluon 
momenta satisfy I gfi. This implies that for I < gfj, the high density 
effective theory becomes non-pcrturbativc and fermion loops in gluon n- 
point functions have to be resummed. The generating functional for hard 
dense loops in gluon n-point functions is given by ^^'^^ 

/AQ pa pP 

where w? = Nfg^iJ.'^/{4:TT^) and the angular integral corresponds to an 
average over the direction of Pq = For momenta I < gfi we have 

to add Lhdl to Lhdet- In order not to over-count diagrams we have 
to remove at the same time all diagrams that become disconnected if all 
soft gluon lines are deleted. Note that the high density effective theory is 
different from the standard hard dense loop approximation, because hard 
loops are resummed only in gluon Green functions, but not in quark Green 
functions or quark-gluon vertex functions. 

The second observation is that the power counting for hard vertices is 
modified by a factor that counts the number of fermion lines in the vertex. 
It is easy to see that four-fermion operators without extra derivatives are 
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leading order (k — 2 — fk = 0), but terms with more than four fcrmion 
fields, or extra derivatives, are suppressed. This result is familiar from the 
effective field theory analysis of theories with short range interactions 

4. Color Superconductivity 

In the last section we saw that hard loops lead to non-perturbative effects 
in the effective theory that require resummation at the scale / ~ gfj,. In ad- 
dition to that, there are logarithmic divergences that have to be resummed 
at exponentially small scales / /iexp(— c/g). The most important effect 
of this type is the BCS instability in the quark-quark scattering amplitude. 
This instability leads to the formation of a gap in the single particle spec- 
trum. We can take this effect into account in the high density effective 
theory by including a tree level gap term 

L = ARf{v- A)tp_ya2T^v + h.c. (9) 

The Dirac matrix F and the angular factor Rf (x) determine the helicity and 
partial wave channel. The magnitude of the gap is determined variationally, 
by requiring the free energy to be stationary order by order in perturbation 
theory. 

At leading order in the high density effective theory the variational 
principle for the gap A gives the Dyson- Schwinger equation 

2g^ f d^q A(g4) n I \ 

where we have restricted ourselves to angular momentum zero and the color 

anti-symmetric [3] channel. is the hard dense loop resummed gluon 

propagator. We also note that equ. (10) only contains collinear exchanges. 
According to the arguments give in Sect. 3 four-fermion operators are of 
leading order in the HDET power counting. However, even though collinear 
exchanges and four-fermion operators have the same power of collinear 
exchanges are enhanced by a logarithm of the small scale. As a consequence, 
we can treat four-fermion operators as a perturbation. 

Sine the electric interaction is screened it is possible to absorb electric 
gluon exchanges into four-fermion operators. At leading order in the high 
density theory the gap equation is completely determined by the collinear 
divergence in the magnetic gluon exchange interaction. This IR divergence 
is independent of the helicity and angular momentum channel. We have 

A(p,) = -|!^ log ( . (11) 
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The leading logarithmic behavior is independent of the ratio of the cutoffs 
and we can set A|| = Aj^ = A. We introduce the dimensionless variables 
variables x = log(2A/((j'4 + e^)) and y = log(2A/(p4 + Cp) where Cq = 
{q1 + A(g4))^/^. In terms of dimensionless variables the gap equation is 
given by 

A{y) = -f^ r dxA{x)K{x,y), (12) 
IStt^ Jo 

where xq = log(2A/Ao) and K{x,y) is the kernel of the integral equation. 
At leading order we can use the approximation K{x,y) = min(a;,y). We 
can perform an additional rescaling x = xqx, y = x^y. Since the leading 
order kernel is homogeneous in x and y we can write the gap equation as 
an eigenvalue equation 

A{y) = xl^ j\xA{x)K{x,y), (13) 

where the gap function is subject to the boundary conditions A(0) = and 
A'(l) =0. This integral equation has the solutions ^® 

A„(5) = A„,osin ( — -^a;o,„x ) , xo,„ = (2n + 1)-^^. (14) 



The physical solution corresponds to n = which gives the largest gap, 
Ao = 2Aexp(— 37r^/(A/2(?)). Solutions with n ^ have smaller gaps and 
are not global minima of the free energy. 

5. Higher Order Corrections to the Gap 

The high density effective field theory enables us to perform a systematic 
expansion of the kernel of the gap equation in powers of the small scale and 
the coupling constant. It is not so obvious, however, how to solve the gap 
equation for more complicated kernels, and how the perturbativc expansion 
of the kernel is related to the expansion of the solution of the gap equation. 

For this purpose it is useful to develop a perturbative method for solv- 
ing the gap equation 20, ^.^^ ^^^w^^ ^j^g kernel of the gap equation as 
K(x,y) = K(){x,y) + 6K{x,y), where Ko(x,y) contains the leading IR di- 
vergence and 6K{x, y) is a perturbation. We expand both the gap function 
A(a;) and the eigenvalue xq order by order 6K, 

A{x) = A(o) (x) + A(i) (x) + A(2) (x) + . . . , (15) 
xo = 4"' + *o'' +x^o^+..., (16) 
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Figure 2. Higher order corrections to the gap equation in the high density effective 
theory. The diagrams shown in this figure correspond to four-fermion operators, fermion 
seff-energy corrections, and vertex corrections. 

where we have defined = g'^XQ/{18TT'^). The expansion coefhcients can 
be found using the fact that the unperturbed solutions given in equ. (14) 
form an orthogonal set of eigenfunctions of Kq. The resulting expressions 
for and A«(x) are very similar to Rayleigh-Schroedinger perturbation 
theory. At first order we have 

4'^ = - ^ (4°))'^'dx^'#A('')(5)5if(a;ox,xoy)A(°)(y), (17) 



_(o) 
J- 

1 - ' 



2k+l 



4'^ = ? 72 dx dyAl^\x)SKixox,xoy)Ak\yl (18) 

Jo Jo 



with A(i)(a;) ^T,Ck^^k\x) and SK = g/{3V2TT)SK. 

We can now study the role of various corrections to the kernel. The 
simplest contribution arises from four-fermion operators. We find 

SK{xox,xoy) = \og{b), b=^^^—j . (19) 

This contribution does not change the shape of the gap function but it gives 
an 0{g) correction to the eigenvalue a;o. This corresponds to a constant pre- 
exponential factor, A*^^) = bA'^*^\ An important advantage of the effective 
field theory method is that this factor is manifestly independent of the 
choice of gauge. The gauge independence of the pre-exponential factor 
is related to the fact that this coefficient is determined by four-fermion 
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operators in the effective theory, and that these operators are determined 
by on-shell matching conditions. 

Another effect that contributes to the eigenvalue at 0{g) is the fermion 
self energy 20,21 ^ one-loop calculation in the high density effective theory 

gives 22,23,24,25,26 

rW^^^.,log(A)^ (20) 

The correction to the kernel of the gap equation is 

- _ _ 9^ - - 

5K{xqx, xoy) = --^ (xox) Ko{xox, xoy), (21) 

and the shift in the eigenvalue is given by 

where (01(5X10) denotes the matrix element of the kernel between unper- 
turbed gap functions, see equ. (17). At this order in g, there is no con- 
tribution from the quark-gluon vertex correction. Note that the quark self 
energy correction makes an 0{g) correction to the eigenvalue, even though 
it is an 0{g^) correction to the kernel. This is related to the logarithmic 
divergence in the self energy. The perturbative expansion of xq is of the 
form 

xo ~ 5log(A) = 0(3°) + Oiglog{g)) + 0{g) + .... (23) 

Brown et al. argued that equ. (22) completes the 0{g) term. At this order 
the spin zero gap in the 2SC phase of TV/ = 2 QCD is 21,27,13 

A = 5127rV.9~^e-^e"^. (24) 

In other spin or flavor channels the relevant four fermion operators are 
different and the pre-exponential factor is modified 5,21,28,29 -^n the CFL 
phase of TV/ = 3 QCD the gap is suppressed by a factor (2/3)^/^2"^/^. 

6. Instanton Correction 

In this section we shall focus on the correction to the perturbatively gener- 
ated gap parameter which is due to instantons. If the density is very large 
instanton effects are exponentially small as compared to perturbative inter- 
actions. In this regime instantons are important for physical observables, 
like the mass of the U{1)a Goldstone mode, that receive no contribution 
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Figure 3. Instanton correction to the perturbatively generated gap in QCD with Nf = 2 
and Nf = 2 + 1 flavors. We show the gap with and without instanton corrections. In 
the Nf = 2 + 1 case we have used m,s = 150 MeV. 



from perturbative interactions, but they make no significant contribution to 
the gap. As the density is lowered instanton effects grow. In the foUowing 
we win use the methods described in the last section in order to determine 
the scale at which instanton contributions to the gap become important. 

Instantons induce a chirality changing four-fermion operator. In QCD 
with Nf = 3 flavors and niu^d <C ms we have 

L = J n{p,fl)dp ^^^^ ^^^ msef,f^eg,g., y-^^{^pR.J^'^pL,g^)(4'R.,f2i^L,g2) 

1 

where the sum over flavors runs over up and down quarks only, /i^2 — gi,2 
(m, d), and the instanton size distribution ji{p, p) is given by 

27V, 



-(V'iiJi (^f^,^1pL,gi)ilpR,f2<^t^u'ipL.,g2] 



{L ^ R) 



(25) 



n{p, p) = C_ 



N 



8^ 
9^ 



exp 



exp 



'Nfp'p^] 



Cn — 



0.466 exp(-1.679iVc)l-34^^ 
(iVc- l)!(Afc-2)! 

= -61og(pA), 



b = — A^c - -Nf . 
3 3-^ 



(26) 
(27) 
(28) 



Note that equ. (25) is an effective interaction for quarks near the Fermi 
surface. In particular, there are no form factors that have to be included. 
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We also observe that the integration over sizes is cut off at p ~ As a 
consequence, instanton effects are of order exp(— 87r^/g^(/i)) ^ {Kqc d / jj)^ ■ 
In addition to the four-fermion vertex given in equ. (25) instantons also in- 
duce a six-fermion operator. This operator has important physical eflFects, 
but it docs not contribute to the gap to leading order in the effective in- 
teraction, so we will not consider it here. In QCD with Nf = 2 massless 
flavors instantons induce a four-fermion operator which can be obtained 
from equ. (25) by making the replacement m^p — > 1. 

We can now compute the instanton correction to the kernel of the gap 
equation. In the three flavor case we find 6Ki = log(6/) with log(6/) = 
9G//g2 and 



V M / V M / V .9^ 

(27r)2 4(7V, + 1) r 



^ ' ^ (29) 



where f3o = llA^c/3 — 2Nf/3 is the first coefficient of the beta function. In 
the two flavor case we get 



IJ'^ \ ^- J \ 9 

(27r)2 4(7V, + 1) 



r(^ 



2(iV2 - 1) TVe 2Nf'''^+^^ ■ 



(30) 



We note that the main difference is a suppression factor (ms/fi) in the 
Nf =3 flavor case. Numerical results are shown in Fig. 3. We observe that 
the instanton correction becomes large for /j, ^ 1250 MeV in Nf = 2 QCD 
and for n ~ 750 MeV in the realistic case of QCD with three flavors and a 
massive strange quark. 

We should note that the fact that the instanton correction to the gap 
is on the order of 100% does not necessarily invalidate the perturbative 
expansion. As explained in the previous section, the quantity that is being 
expanded is not the gap A, but log(/z/A). In Fig. 4 we compare the size of 
the leading order 0{l/g) term, the sub-leading 0(log(g)) and 0(1) terms, 
as well as the instanton contribution to the logarithm of the gap in Nf = 3 
QCD. We observe that exponentially small instanton contributions start to 
dominate over all other contact terms for jjL ~ 600 MeV. 
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Figure 4. This figure shows successive terms in the perturbative expansion of log(/^/A) 
in QCD with Nf = 3 flavors. 

For ij, ~ 500 MeV the instanton term becomes comparable to the leading 
0{l/g) term and we can no longer consider instantons effects to be a small 
correction. In this regime it probably makes more sense to do an instanton 
calculation and consider perturbative gluon exchanges as a correction. This 
is the approach originally suggested in We should note, however, that 
for n < 500 MeV the instanton calculation requires some phenomenological 
input because the instanton size distribution equ. (26) is no longer reliably 
calculable. 

Ideally, we would like to check perturbative calculations of the gap 
against numerical calculations on the lattice. While this cannot be done in 
the realistic case of QCD with Nc — 3 colors, the comparison is possible for 
QCD with Nc — 2 colors and an even number of flavors, or for QCD with 
Nc = 3 colors and a non-zero chemical potential for isospin rather than 
baryon number. In both of these cases it is also possible to separate the 
perturbative and instanton contributions by measuring both the gap and 
the mass of the U{1)a Goldstone boson 

7. Numerical Estimates 

In this section we shall address a few more issues related to the magnitude of 
the gap. We will estimate the size of certain higher order corrections using 
numerical solutions of the gap equation. The first problem we wish to study 
is the importance of the fermion self energy correction. In Sect. 5 we saw 
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Figure 5. Numerical solution of the perturbativc gap equation. We show the results 
with and without the fermion self energy included. We also show the asymptotic result, 
which corresponds to a reduction of the gap by a factor exp(— (tt'^ + 4)/8) ~ 0.18. 



that at asymptotically large chemical potential non-Fermi liquid effects in 
the fermion self energy reduce the gap by a factor exp(— (tt^ -|-4)/8) ~ 0.18. 
In order to study the problem at moderate densities we consider the gap 
equation 



A(p4 



dq4 Z{q4) 
\/qi 



A(94) 



A(g4)- 



log 



with 



^(94) 



1 



97r- 



■log 



94 



(31) 



(32) 



and 6o — 2567r^(7~^. If the density is very large then the scale inside the 
logarithm in equ. (32) does not matter. For our numerical estimates we have 
used the electric screening mass which is the scale suggested by one-loop 
calculations 

Numerical results are shown in Fig. 5. We observe that the gap is indeed 
reduced by the effects of fermion wave function renormalization, but at 
moderate density the reduction is smaller as compared to the asymptotic 
result. This is related to the fact that for ^ ^ 500 MeV and A 50 MeV 
the logarithm log(m£)/A), which is formally ©(l/g), is numerically not 
very large. 
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Figure 6. Numerical solution of the gap equation with running coupling constant effects 
included. We show the result with a fixed coupling g = g{li), with a coupling that runs as 
a function of the momentum transfer as long as the momentum is larger than the electric 
screening scale, and a coupling that runs as a function of the momentum without any 
restrictions. 



The second issue we wish to address is the role of running couphng con- 
stant effects in the gap equation . For the numerical estimates presented 
in the last section we have used the weak coupling result equ. (24) with the 
coupling constant evaluated at the scale /i. It is easy to see that variation 
in the scale correspond to 0{g^) corrections, which is of the same order as 
other terms that have been neglected. 

In order to estimate the size of these effects we consider the gap equation 

n 4/(3/3o 



i 187r2 



M 



A3 



2 log 



(33) 



where Am = (7rm|,|p4 ± q4|/4)^/'^ is the scale that characterizes magnetic 
gluon exchanges. For momenta above the electric screening scale we have 
used the one-loop running coupling at the scale set by the momentum 
transfer. For momenta below the screening scale the coupling is frozen. In 
the language if the high density effective theory this means that the we have 
performed the matching at the electric screening scale. The four-fermion 
operator acquires an anomalous dimension which is equal to the QCD beta 
function . Below the electric screening scale the gluonic interaction is 
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effectively three-dimensional and does not run. 

In Fig. 6 we show results obtained by solving equ. (33) numerically. 
For comparison we also show results obtained by solving a gap equation 
in which the coupling is allowed to run all the way down to the magnetic 
scale {rnjjAy^^ . This approximation was suggested by Beane et al. . We 
observe that for a moderate chemical potential ~ 500 MeV the running 
coupling constant does not lead to very large effects. The reason is that 
there is no large hierarchy between the scales mn and /x. If the chemical 
potential is very large, /i ^ 10 GcV, the gap is increased by about 50%. 
This effect slowly disappears at asymptotically large chemical potential. 
The situation is different in case of the gap equation proposed by Beane et 
al. In that case the gap equation involves an extra large logarithm log(/Lt/A) 
and the pre-exponential factor in the asymptotic solution is modified 

8. Conclusions 

We discussed an effective field theory for QCD at high baryon density. We 
studied, in particular, the problem of power counting in the high density 
effective theory. We showed that the power counting is complicated by 
"hard dense loops" , i.e. loop diagrams that involve the large scale and 
proposed a power counting that takes these effects into account. The mod- 
ified I counting implies that hard dense loops in gluon n-point functions 
have to be resummed below the scale gn, and that four fermion operators 
are leading order in the HDET power counting. 

Wc used the high density effective theory to study the size of instanton 
corrections to the gap in superfluid quark matter. We found that instanton 
effects are very large in the regime /i ~ 500 MeV which is of physical 
interest. We argued that numerical calculation in QCD with Nc = 2 colors, 
or QCD with = 3 colors and non-zero isospin chemical potential, will 
of great help in determining the gap in TVc = 3 QCD at non-zero baryon 
density. We also studied a renormalization group improved gap equation. 
We find no significant corrections as compared to a gap equation with a 
fixed coupling. 

Acknowledgments: This work was supported in part by US DOE grant 
DE-FG-88ER40388. 
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